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Abstract 

This paper gives a new integrable hierarchy of nonlinear evolution equa- 
tions. The DP equation: nit + unix + Smux = 0, m = u — Uxx, proposed re- 
cently by Desgaperis and Procesi 0, is the first one in the negative hierarchy 

2 2 

while the first one in the positive hierarchy is: rrit = 4:{m~'s)x — 5{m~3)xxx + 

2 

{m^s^xxxxx- The whole hierarchy is shown Lax-integrable through solving a 
key matrix equation. To obtain the parametric solutions for the whole hier- 
archy, we separatedly discuss the negative and the positive hierarchies. For 
the negative hierarchy, its 3 x 3 Lax pairs and corresponding adjoint repre- 
sentations are nonlinearized to be Liouville-integrable Hamiltonian canonical 
systems under the so-called Dirac-Poisson bracket defined on a symplectic 
submanifold of M^^. Based on the integrability of those finite-dimensional 
canonical Hamiltonian systems we give the parametric solutions of the all 
equations in the negative hierarchy. In particular, we obtain the parametric 
solution of the DP equation. Moreover, for the positive hierarchy, we consider 
the different constraint and use a similar procudure to the negative case to ob- 
tain the parametric solutions of the positive hierarchy. In particular, we give 

2 2 

the parametric solution of the 5th-order PDE mt = 4(m~3)^ — 5{m~3)xxx + 

2 

{m~3)xxxxx- Finally, we discuss the stationary solutions of the 5th-order PDE, 
and particularly give its four peaked stationary solutions. The stationary so- 
lutions may be included in the parametric solution, but the peaked stationary 
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solutions not. The 5th-order PDE does not have the cusp sohton although it 
looks like a higher order Harry-Dym equation. 

Keywords Hamiltonian system, Matrix equation, Zero curvature representation, 
Integrable equation. Parametric solution. Peaked stationary solution. 
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1 Introduction 

The inverse scattering transformation (1ST) method is a powerful tool to solve inte- 
grable nonlinear evolution equations (NLEEs) . This method has been successfully 
applied to give soliton solutions of the integrable NLEEs. These examples include 



the well-known KdV equation [16|, which is related to a 2nd order operator (i.e. Hill 
operator) spectral problem [|17[ [19|, the remarkable AKNS equations JI], Q], which is 
associated with the Zakharov-Shabat (ZS) spectral problem [Q, and other higher 
dimensional integrable equations. 

To look for as many integrable systems as possible has been an important topic 
In the theory of integrable system. Kaup |jl3|] studied the inverse scattering problem 



for cubic eigenvalue equations of the form i/j^xx + GQi^x + 6Rip = Xip, and showed a 



5th order partial differential equation (PDE) Qt + Qxxxxx + 30{QxxxQ + \QxxQx) + 
ISOQxQ^ = (called the KK equation) integrable. Afterwards, Kuperschmidt [|l^ 
constructed a super-KdV equation and presented the integrability of the equation 
through giving bi-Hamiltonian property and Lax form. In 1984 Konopelchenko and 
Dubrovsky [Q presented a 5th order equation ut = {u~'^^^)xxxxx and pointed out 
that this equation is a reduction of some 2-1-1 dimensional equation. We have found 



the parametric solution and some traveling wave solutions of this equation ||TT . 

Recently, Degasperis and Procesi proposed a new integrable equation: rrit + 
urux + ^mux = 0, m = u — Uxx, called the DP equation, which has the peaked soliton 
solution. The DP equation is an extention of the Camassa-Holm (CH) equation 
[|], and is proven to be associated with a 3rd order spectral problem |0: ipxxx = 
ipx — Xmi/j and to have some relationship to a canonical Hamiltonian system under 
a new nonlinear Poisson bracket (called Peakon Bracket) 

The purpose of the present paper has two folds: 

• extend the DP equation to a new integrable hierarchy of NLEEs through 
studying the functional gradient of the spectral problem and a pair of Lenard's 
operators; 
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• connect the DP equation to some finite-dimensional integrable system, give its 
parametric solution from the view point of constraint, and furthermore study 
the parametric solution, stationary solution and cusp soliton solution. 

The whole paper is organized as follows. Next section is saying how to connect 
a spectral problem to the DP equation and how to cast it into a new hierarchy of 
NLEEs, and is also giving the pair of Lenards operators for the whole hierarchy. 
In section 3, we construct the zero curvature representations for this new hierarchy 
through solving a key 3x3 matrix equation, and therefore this hierarchy is inte- 
grable. In particular, we give the matrix Lax pair of the DP equation, which is 
equivalent to the form in Ref. 0, as well as the Lax pair for a 5th-order equation 

2 2 2 

rrit = 4(m"3)^ — 5{m~^)xxx + {itT'~'^)xxxxx- We will see that the DP equation is 
included in the negative hierarchy while the 5th-order equation in the positive hi- 
erarchy. To obtain the parametric solutions for the whole hierarchy, we separatedly 
discuss the negative and the positive hierarchies. In section 4, we deal with the neg- 
ative hierarchy. Its 3 x 3 Lax pairs and corresponding adjoint representations are 
nonlinearized to be Liouville-integrable Hamiltonian canonical systems under the 
so-called Dirac-Poisson bracket defined on a symplectic submanifold of . Based 
on the integrability of those finite-dimensional canonical Hamiltonian systems we 
give the parametric solutions of the all equations in the negative hierarchy. In par- 
ticular, we obtain the parametric solution of the DP equation. Section 5 copes with 
the positive hierarchy. We consider the different constraint between the potential 
and the eigenfunctions. Under this constraint the 3x3 Lax pairs and corresponding 
adjoint representations of the positive hierarchy are nonlinearized to be Liouville- 
integrable Hamiltonian canonical systems in the whole M^^. Then we obtain the 
parametric solutions of the positive hierarchy. In particular, we give the parametric 

2 2 2 

solution of the 5th-order PDE nit = 4(m~3)^ — 5{m~3)^^^ -|- (m~3)^^^^^. Finally, in 
section 6 we discuss the stationary solutions of the 5th-order PDE, and particularly 
give its four peaked stationary solutions. The stationary solutions may be included 
in the parametric solution, but the peaked stationary solutions not. The 5th-order 
PDE does not have the cusp soliton, either. 

2 Spectral problems and Lenards operators 

Let us consider the following spectral problem 




a 
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and its adjoint problem 

r... = + AmV^*, (2) 

where A is a spectral parameter, m is a scalar potential function, ■0 and tp* are the 
spectral wave functions corresponding to the same A, ip* is not conjugate of ip, and 
a = constant. Then, we have 

where 

poo 

E = — mip'ilj*dx = const. 



Here during our computation about the functional gradient ^ of the spectral pa- 
rameter A with respect to the potential m, we need the boundary conditions of 
decaying at infinities or periodicity condition for the potential function m. A gen- 
eral calculated method can be seen in Refs. ||^, ^ . 
Now, we denote the function ViA by 

V.A=^'»^--^'^-\ (4) 
Then, we have the following equality: 

^(VA, ViA)^ = A J(VA, ViA)^ 
where K and J are two matrix operators 



K 



J 



Ad + 5a'^d^ - a^d^ 

-^d 



?)a^{2md + dm) 

md + 2dm 



d = — 

dx ' 

or we rewrite Eq. (5) as the following Lenard spectral problem form 

KVX = AVVA, (5) 

where 

K = 4d-5a^d^ + a^d^, 

J = 3a%2md + dm){-a^d^ + d)-\md + 2dm). 



4 



Without loss of generality, we assume a = 1 below. Then, K, J read 

K = Ad-5d^ + d^, (6) 
J = 3{2md + dm){d-d^)~\md + 2dm). (7) 

Hint: Here we do not care about the Hamiltonian properties of the operators 
K, J, but need 



27 ^ ' 



They yield 



jC-^ = R-^J = 3{d - d^y\4: - d'^y\2md + dm){d - d^y\md + 2dm), 
which are actually the two recursive operators we need in the next section. 



3 Zero curvature representations and an integrable 
hierarchy 

Letting ijj = i/ji, we change Eq. (|l]) to a 3 x 3 matrix spectral problem 

= f/(m,A)*, (8) 

(9) 

Apparently, the Gateaux derivative matrix f/^, {^) of the spectral matrix U in the 
direction ^ G C°°(M) at point m is 





( 





1 


o\ 






U{m,X) = 










1 








V 


—mX 


1 


J 




J 



de 



e=0 



/ \ 


\ -CA / 



(10) 



which is obviously an injective homomorphism. 

For any given C°°-function G, we construct the following 3x3 matrix equation 
with respect to = V{G) 



- [[/, V] = U,{KG - AVG). 
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Theorem 1 For the spectral problem ^ and an arbitrary C°° -function G, the ma- 
trix equation ^Tl\ ) has the following solution 

/ rc + sAae-iTG sGa^-sxe-^rc ~6G \ 

V = \\ rGa: + 3X{d^e-^TG + 2mG) ~2{G ~ G^^x) -SG^ - SXQ-^TG 

\ TGxx + iX{d + Xrn)e-^TG -QG - 3X{d-^TG - 2mG) -2G - G^x - iXdQ-^TG j 

(12) 

where 9 = T = md+2dm, T = A-d'^. Therefore, J = 3T*e-^T, K = TO. 



Proof: Let 

/ Vu 1^13 \ 

V= V21 V22 V23 , 

\ ^31 V32 Vs3 J 

and substitute it into Eq. ([TT|). That is a over-determined equation. Using some 
calculation technique |2^, we obtain the following results: 

V12 = 3XG^-3X^e-^TG, 
Vi3 = -6XG, 

V21 = XrG^ + 3X\d^Q-^TG + 2mG), 

V22 = —2X{G — Gxx), 

V23 = -3XGx-3X^Q-^TG, 

V31 = XrGxx + 3X\d + Xm)Q-^TG, 

V32 = -XQG-3X\d-^TG-2mG), 

V33 = -2XG - XGxx-3X^dQ-^TG, 



which completes the proof. 

Theorem 2 Let Go E Ker J = {G e G^ 



JG = 0} and G_i G Ker K 



{G e 



Then, 



KG = 0}. We define the Lenard's sequences as follows 



G, = C^Go = C^^^G-i, J e Z. 



(13) 



1. the all vector fields Xk = JGk, A; G Z satisfy the following commutator repre- 
sentation 



Vk,x-[U,Vk] = U,{Xk), VA;gZ; 



(14) 
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2. the following hierarchy of nonlinear evolution equations 

mt, =Xk = JGk, yk e Z, (15) 
possesses the zero curvature representation 

Ut,-Vk,^ + [U,Vk]=0, VA;gZ, (16) 

where 



> 0, 

. -E7i., k<o, 

and V{Gj) is given by Eq. ( |7^ with G = Gj. 
Proof: 

1. For A; = 0, it is obvious. For A; < 0, we have 

-1 

Vk,.-[U,Vk] = -Y,{V,{G,)-[U,ViG,)])X'^'-^-'^ 

j=k 
-1 

= -Y^U, [KG, - X'KG,^,) A2('=-^-i) 

j=k 

\j=k 

= U^KGk-i-KG^.X^') 
= U,{KGk-i) 
= U^Xk). 

For the case of A; > 0, it is similar to prove. 

2. Noticing Ut,^ = U^:{mt^), we obtain 

Ut, - Vk,. + [[/, Vu] = U,{mt^ - Xk). 
The injectiveness of U^, imphes item 2 holds. 



So, the hierarchy (p!5| ) has Lax pair and is therefore integrable. In particular, 
through choosing G_i = — | G Ker K, ([T5|) reads 

= -J£'=+i ■ ^, A; = -1,-2,..., (18) 



where C = J ^K. Set m = u — u^x, then it is easy to see the first equation in the 
hierarchy is exactly the DP equation M 



rrit + urrix + 3mUx = 0, t = t_i. 

This equation has the following Lax pair: 

^x = t/(m,A)^, 
= V{m,X)^, 

where U{m, A) is defined by Eq. (^), and V{m, A) is given by 

/ w^ + |A-^ -u -A-i \ 



(19) 



(20) 



U 



V{m, A) 

which can be changed to the form in Ref. 



—u 



—U 



\y Ux + umX —Ux — |A ^ 



(21) 



V't + A Vxa' + - ( + 2 A ^ ) = 0, = ^Z-i. 



(22) 



2 _ 

Let us choose a kernel element Go from KerJ: Go = mT^. Then Eq. (^) reads 
the following integrable hierarchy 



nit^^ = JC^m k = 0, 1,2, .... 
In particular, the equation 

nit = 4(m"5)a; - b^m'^xxx + {m~^)xxxxx, 

has the Lax pair: 

^x = f/(m,A)^, 
where U{m, A) is defined by Eq. (^), and Vi{m, A) is given by 



(23) 



(24) 



(25) 



Tm 3 3(to 3 )x 



— 6to 3 
-3(m~i)j; 
-2m 3 — (to 3 jj.^ 



(26) 



with the operator F = 4 — 9^. Eq. (2^) is therefore a new integrable equation. 

In the next two sections we will give parametric solutions for the negative order 
hierarchy (|18D and the positive order hierarchy (^). 
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4 Parametric solution of the negative order hier- 
archy ( [tBl ) 

To get the parametric solution, we use the constrained method which leads finite 
dimensional integrable systems to the PDEs. Becasue Eq. is a 3rd order eigen- 
value problem, we have to investigate itself together with its adjoint problem when 
we adopt the nonlinearized procedure ||^. Ma and Strampp |18[ ever studied the 
AKNS and its its adjoint problem, a 2 x 2 case, by using the so-called symmetry 
constraint method. Now, we are discussing a 3 x 3 problem related to the hierarchy 
(|T5]). Let us return to the spectral problem and consider its adjoint problem 



mX\ 




( \ 


-1 -1 




^2 


-1 y 




\rs J 



(27) 



where i/j* = ip^^ 



3- 



4.1 Nonlinearized spectral problems on a symplectic sub- 
manifold 

Let Xj (j = 1, A^) be N distinct spectral values of (|) and (p7|), and qij,q2j,Q3j 
and Pij,P2j,P3j be the corresponding spectral functions, respectively. Then we have 

qix = 92, 

q2x = gs, (28) 
qsx = -mkqi + g2, 

and 

Pi^ = mAp3, 

P2x = -Pi-P3, (29) 
Psx = -P2, 

where A = diag{Xi, \n), qk = {qki,qk2, ■■■,qkNY , Pk = {Pki,Pk2, ■■■,PkNV , k = 
1,2,3. 

Now, we consider the following (6A^ — 2)-dimensional symplectic submanifold in 

M = {{p,qf eR^^l F = 0, G = 0} (30) 

where p = (pi,P2,P3)^, q = (gi,g2,?3)^, F = {Aqi,ps)-1, G = (Ag2,P3)-(Ag3,P2) + 
{Aq2,pi), and (■, ■) stands for the standard inner product in M^. 
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When we restrict the above two systems in to the submanifold M, we obtain 
a constraint of m relationship to the spectral function p, q: 

^_o (Ag2,P2)-(Ag3,Pi+P3) .3^^ 

(A2g2,P3) + (A'gi,P2) 

Notice: {^^q2,'Pz) + (^^^1,^2) 7^ is necessary because it assures M is a sym- 
plectic submanifold in M^^. 

Under the constraint (^) the two systems (|28|) and (129|) are nonlinearized as 
follows: 

qix = q2, 

q2x = gs, (32) 



„ _ O (Ag2,P2)-(Ag3,P3)-(Ag3,Pl) A ^ _|_ « ■ 



and 



„ _ r) (Ag2,P2)-(Ag3,p3)-(Ag3,pi) a 

y^lx - ^ (A2g2,p3) + (A2qi,p2> ^^^3' 

P2x = -Pl-P3, (33) 
PSx = -P2- 

They are forming a (6A^ — 2) dimensional nonlinear system on M with respect 
to p, q. Is it integrable? To see this, in we modify the usual standard Poisson 
bracket of two functions Fi, F2 as follows: 

which is still antisymmetric, bilinear and satisfies the Jacobi identity. 
Obviously, 

{F,G} = (A'g2,P3> + (A2gi,p2> ^ 0. (35) 

Because we are discussing the system on the submanifold M, we need to intro- 
duce the so-called Dirac- Poisson bracket of two functions /, g on M: 

{/, g}D = {/, 9} + ({/, F}{G, g} - {/, G}{F, g}) (36) 

which is satisfying the Jacobi identity. 
Now, we choose a simple Hamiltonian 

H = {q2,Pi+P3) + {q3,P2) , (37) 
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then, the two systems (^) and (|33D have the canonical Hamiltonian form on M: 



Qij,x = {qij, h}d, 
Q3j,x = {q3j, H}d; 

Plj,x = {pij,H}D, 
P2j,x = {p2j,H}D, 
P3j,x = {p3j,H}D, 

which can be in a brief form rewritten as: 

Qx = {q,H}D, 
Px = {p, H}d. 

In this calculation procedure, we have used 

{H,G} = 2(^{Aq2,P2)-{Aq3,Pi+P3)), 
{H,F} = 0. 

It is easy to check that = 0, i.e. H is invariant along the flow 
H = Cd {Cd is constant) along this flow and 

M = ^ ((gi,P2) + (g2,P3)) - ]^H, 

then we have 



(3J 



(39) 



Assume 



(40) 



u — 



(41) 



which is exactly related to the DP equation (|I^). 

To show the integrability of canonical system (0), we need to consider the 
nonlinearization of the time part of the Lax representations. 

4.2 Nonlinearized time part on this submanifold 



Let us turn to the time part (|20| ) of the Lax pair for the TD equation (|19|). Then 
the corresponding adjoint problem reads 





^ —Ux — |A ^ —u 


—Ux — umX ^ 






( \ 




u \\-^ 







^)* = 






^ u 


Ux + \\-^ ) 






Us / 



(42) 



11 



We also consider the constrained system of the time part on M. Thus, under the 
constraints (|31|), (|40|), and 

Ux = ^ ((gs - gi,P3) - , (43) 



Eqs. (|20D and (^) are nonhnearized as: 



qit = u^qi + 3A ^qi - uq2 - A ^ga, 



q2t = uqi - I A ^q2 - uq^, (44) 
qzt = u^qi + umKqi - \k-^qz - u^qs] 



and 



Pit = -UxPi - fA^Vi - up2 - u^Ps - umAp3, 

P2t = upi + |A"V2, (45) 

P3t = A^Vi + up2 + UxPs + lA^Vs, 

respectively, where each qk,Pk and A are the same as section 2, and A~^ is the inverse 
of A. 
Let 

2 11 

I = -{A-'qi,pi)--{A-'q2,P2)--{A''q3,P3)-{A''q3,Pi) 

+ \{quP2f-\{q2,P3f-lH{q,,p2) 

-u{{q2,Pi) + {q3,P2)-H)-ul (46) 

then the two systems (Q) and ( ^5] ) are expressed in a canonical Hamiltonian form 
on M: 



qt = {qJ}D, 
Pt = {p, I}d, 

where p = {puP2,P3V,q = (<?!, ^2, ^s)^- 

In the above calculations, we use the following equalities: 

F = (Agi,p3)-1 = 0, 

G = (Ag2,Pi +P3) - (Ag3,P2) = 0, 

= {Aq2,P3) - {Aqup2) = 0, 

F^^ = {Aqs,p3) -2{Aq2,P2) + {Aqupi) + 1 = 0, 

{G,I} = 2u(^{Aq2,P2) - {Aq3,pi+P3) 
{F,I} = 0. 



(47) 



12 



Theorem 3 



{i/,/}^ = 0, (48) 
that is, two Hamiltonian flows commute on M . 
Proof: By the definition, 

{H., = {H, 1} + {{H, F}{G, 1} - {H, G}{F, I}) , (49) 

we need to compute each Poisson bracket in this equality. 

{H,l} = Mx((P2,?i -?3) - (?2,Pi)) + M((P3,gi -ga) - (gi,pi)) = 0, 
{H,F} = (Agi,p2)-(A?2,P3) = (-(Agi,P3) + l). = 0, 
{F,I} = 2u,{{Aqi,p3) -l) + u{-{Aqi,p3) + l). = 0, 

complete the proof. 

By Theorem 2, we know that the hierarchy (^) has the Lax representation: 

= ?7(m,A)*, (50) 

-1 

= -5^V(G'fc)A2(^-^-i)vi>, ^■<0, (51) 

k=j 

where V{Gk) is given by Eq. (0) with G = Gk- 

In last section, we have investigated the nonlinearized systems of spectral prob- 
lem and its adjoint. Now, we study the nonlinearizations of time part (^) and its 
adjoint problem: 

-1 

k=j 

where V'^{Gk) is the transpose of V{Gk)- 

Let = (gifc,g2fc,g3fc)^, "^l = {Pik,P2k,P3kV be the eigenfunctions correspond- 
ing to the eigenvalues {k = 1,...,N) of (§) and (|27|) . Let us start from the 



constraint G_i = —^Ylf=i^^j- This constraint is giving the symplectic subman- 
ifold M we need. Let the two antisymmetric operators act on this constraint, we 
have: 

G, = -^{A'^+\,Ps), 3<0. (53) 
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Therefore, a comphcated calculation yields the following formulations: 

1 
6 



TG 



3,^ 



TG 



J,XX 



QG. 



d-'TGj - 2mGj 



d'^Q-^TGi + 2mGi 



1{{A'^''%,P2) + {A'^'-'q2,Pi) - {A'^''%,P2)), 
l[m[{A'^^%,p,) + {A'^^%,Ps)) 
+2{A'^^\,,p,) + {A'^^'q,.p,) - (A^^+V,Pi +P3> 
\{{A'^^'q2.P.+P,) - {A''^\,.p,)), 
-\{{A''^'q2,Pi+P,) + {A''^%.P2)). 
\{{A'=^'q2,Pi) + {A'^'-'q..P2) - (A^^+V,P2>), 



Substituting these equalities to Eqs. (|5lD and (|5^ , with a similar computational 
method to section 3 we find the nonlinearized systems of the time part ( |5T[ ) and the 
adjoint time part (|52|) are cast in a canonical Hamiltonian system on the (6A^ — 2)- 
dimensional symplectic submanifold M: 



Ptk = {P^Fk}D, 



k 



(54) 



where 



Ql,Pl , 



+ E [ - 1^ (A^^+^gi,Pi) (A2('=-^)-igi,pi) - i (a2^+3<?2,p2) (A2('=-^)-ig2,p2) - i (A2^+^g3,P3) (A2('=-^)-ig3,P3) 

+ 1 (A2^+3q2,P2) (A^C^-^)-!,!,^!) - 1 (A2^+3g3,p3) ( A^C^"^) "l^i , pi ) + 1 (A^^+S 53,^3) (A^C^-^)-!^^,^^) 
-i (A2^+^gi,P2) (A2('=-^)-lgi,P2) + i (A^'+'quPs) (A^^'-'^-'quPs) + \ {A^'+\2,P3) (A2(*-^)-lg2,P3) 
+ i ( (A2^+3gi,p2> - (A2j+3g2,P3) ) ( (A2('=-^)-lg3,P2) - (A2('=-^)-lq2,pi) ) 
+ i (A2^+^qi,P3) ( {f^'^^-^^-\„p,) - 2 (A2('=-i)-lg3,pi) - {a^(^-^)-\^,p^)^ 

+ E [- (A2(*-^)gi,pi) - i(A2^+2q3,p3) (A2('=-^)g3,p3) 

j=fe 

+ 1 (A2^+2gi,p2) (A^C'-^jgi.pa) - 1 (A2:'+2gi,p3) ( A2{*-^) , P3) - 1 ( A^^+^g^ , P3) ( A^C^"^) 53 , P3) 

-i((A2^+2gi,P2)+(A2^+2g2,P3))((A2('=-^)g3,P2)+(A2('=-i)g2,Pi)) 

+ i {A^'+\UP-S) (A2('=-^'q3,P3) + i (A2^+2g3,P3) (A2{^-^)gi,pi) - i (A2^+2gi,p3) ( A2{^-^) 51 , p^ ) 
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Apparently, F_i = I. Furthermore, by the Dirac-Poisson bracket (|36|) on sub- 
manifold M we obtain the following theorem. 



Theorem 4 All Hamiltonian flows ( \3^ ) and ( \54) commute on M. 
Proof: Through a lengthy calculation, we have 

{H,Fk}D = 0, 

which completes the proof. 

Therefore, all Hamiltonian flows (|54D are integrable on M. Particularly, the 
Hamiltonian system (|39D is integrable. 



Remark 1 In the calculation procedure and the proof of Theorem^ and Eq. (^), 
we used the following facts: 

{'?i,i^} = te,i^} = 0, 

{Pi,i^} = {P2,F} = 0, 
{p^,F} = -Ap3, 
{F,Fk} = 0, 

{G,F,} = ((Aq2,P2) - {Aq3,pi+pM{^""'\,P2) + {A''^'q2,P3)), 
and on M all of these equalities F = G = F^ = F^x = are valid. 

4.3 Parametric solution 

Theorem 5 Let p{x,tk),q{x,tk) (p{x,tk) = iPi,P2,P3V , qix,tk) = iqi,q2,q3V, k = 
— 1,-2,...) be the common solution of the two integrable flows ^3^) and ^5^), then 



_ ^ (Ag2,P2)-(Ag3,Pi+P3) 

satisfies the negative order hierarchy ([7^. 
Proof: Noticing the following formulas 

Gk = -^(A^'^+V.Pa), ^ = -1,-2,-3,..., 
T = md + 2dm, T* = dm + 2md, Q = d — d^, 
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and Eq. (|5^) , we directly compute and find Eq. (^) satisfies rrit^ = 3T*6 ^TGk 
which completes the proof. 

In particular, we obtain the following theorem. 

Theorem 6 Let p{x,t),q{x,t) (p{x,t) = {pi,p2,pzf ,q{x,t) = {qi,q2,q^f) he the 
common solution of the two integrable flows and (j^^, then 

^ ^ r, {^<l2{x,t),p2{x,t)) - {Aq3{x,t),pi{x,t) +p3{x,t)) 

{A^q2{x,t),ps{x,t)) + {A^q,{x,t),p2{x,t)) ' ^ ^ 

1 1 

u = -i{qi{x,t),p2ix,t)) + {q2{x,t),ps{x,t)))--H, (57) 

satisfy the DP equation: 

rrit + urrix + "imUx = 0. (58) 



Proof: Let 



A = (Ag2,P2) - (Ag3,Pi +P3) , 

B = {A^q2,P3) + {A^qi,P2), 

C = {A'^q3,P3) - (A^i^pi+ps) . 



Then through a lengthy calculation we have 



At = uG + uraC — 2uxA = umC — 2uxA, 

Bt = G-uC + u^B = -uC + u^B, 

Ax = -2G -mG = -mC, 

Bx = G. 

By the above equalities, we obtain 

rrit + unix + 3mUx 

= ^ [AtB - ABt + u{AxB - ABx) + ?>UxAB 
2 / 

= — [umBG + uAG + uAxB - uABx 

= '^(B{mG + Ax) + A{G - Bx)) 
= 0. 

In the above proof procedures, we have used the following equalities: F = G 
0, Fx = Fxx = on M. 
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Similarly, we can discuss the parametric solution of the positive order hierarchy 
(PSD- That needs us to consider a new kind of constraint and related integrable 
system, which we deal with in the next section. 



5 Parametric solution of the positive order hier- 
archy ( [23l ) 

Let us directly consider the following constraint: 

N 

Go = Y,E,VX„ (59) 
i=i 

where Ej'VXj = Xjqijp^j is the functional gradient of Xj for the spectral problems 
(H) and (p7|), and qkj, Pkj {k = 1,2,3) are the related eigenfunctions of Xj. Then 
Eq. (^) is saying 

m={Aqi,p3)~^ (60) 

which composes a new constraint in the whole space R^^. Under this constraint, the 
spectral problem and its adjoint problem (|27| ) are able to be cast in a Hamiltonian 
canonical form in M^^: 

(H^): '^ = 1''^?' (61) 
Px = Ip, H+}, 

with the Hamiltonian 

H+ = {q2,pi+P3) + {q3,P2)+ .^'^ . (62) 

V(Agi,P3) 



To see the integrability of the system (|6TD, we take into account of the time part 
= Vi{m, A)\E' and its adjoint = —V^{m, A)\E', where Vi{m, X) is defined by Eq. 
(pGl). Under the constraint (0), the time part and its adjoint are also nonlinearized 
as a canonical Hamiltonian system in E'" 

(F,): "'"i"'!!'!' (63) 

Pt = {p,Fi}, 

with the Hamiltonian 
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(A'gi,P2) + (A'g2,P3) 
V(Agi,P3) 

-^((Agi,Pi)' + 4(A(/2,P2)'+ (Ag3,P3)') 

+ ^((Agi,p3)' + (Ag2,P3)' - (A9i,P2)') 

+3 (Agi,P3) ((A9i,Pi) - 2 (Ag3,Pi) - (Ag3,P3)) 

+2 (Agi,pi) (Ag2,P2) + 2 (Ag2,P2) (Ag3,P3) - i^QuPi) (Ag3,P3) 

+3((Agi,p2) - (Ag2,P3))((Ag3,P2) - (Ag2,Pi)). (64) 



After a calculation of the Poisson bracket of {H^,Fi}, we know that the two 
canonical Hamiltonian flows commute, i.e. 

{H^,F,} = 0. 

Furthermore, under the constraint (0) the nonlinearized systems of the general time 
part = ^JZg y(Gj)A2(^--'-i)^ (A; > 0, k e Z) and its adjoint problem produce 
the following canonical Hamiltonian system in 



= {Q,Fk}, 



k — 1,2, 3, 



(65) 



with the Hamiltonian 



^ (A^Nl,P2)+(A^N2,P3) 
3 = 

+2(A2^+ig2,P2)(A2{'=-^)-igi,pi)-(A2^+ig3,P3)(A2('=-^)-i<?i,pi)+2(A2^+i<?3,P3)(A2(''-^)-i<?2,P2) 

(A2^+igi,P2) (a^('=-^'-^<?i,P2) + I (A2^+i<7i,P3) (a2('=~^)-i<?i,P3) + ^ (A2^+^2,P3) (A2(*-^)-ig2,P3) 
+3((A2^+igi,P2) - (A2^+ig2,P3) ) ( (a2('=-^)-153,P2) - (A2('=-^)-ig2,Pi) ) 

+3 (A2^ + lgi,p3) ( (A2{'=-^)-lgi,pi) - 2 (A2('=-i)-l53,pi) - (A2('=-:')-lg3,p3) j 

+ E [ - i (A'^9i,Pi) (A2(*-^)gi,pi) ^ ^ (A2^g3,P3) (A2('=-^)g3,P3) 

+ ^ (A2^gi,P2) (A2('=-i)gi,p2) - ^ (A2^gi,P3) (a2('=~^)<7i,P3) - ^ (A^:- 92,^3) (a2('=-:')<72,P3) 
-3(^ (A2^gi,P2) + (A2^g2,P3) ) (A2('=-^)g3,P2) + (a2('=-^)52,Pi) j 

+3(A2jgi,P3)(A2(fe-i)g3,p3^ +3(A2jq3,P3) ( A^C^^^) qi , pi ) - 3 ( A^J gi , p3 ) (| A^C^"^) gi , pi^ 
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Apparently, when k = 1, is exactly Eq. (|64D. Furthermore, through a lengthy 
computation, we obtain 



{i/+,Ffc} = 0, A; = 1,2,..., (66) 

which represents each Hamiltonian t-fiow {Fk) commutes with Hamiltonian x-fiow 
{H'^). Thus, all Hamiltonian canonical systems (Fk) are integrable in M.^^ . Partic- 
ularly, the nonlinearized spectral problems is integrable. 
Like last section, we also have a similar theorem 

Theorem 7 Let p{x,tk),q{x,tk) (p{x,tk) = {pi,p2,P3)'^ , q{x,tk) = ^2, gs)'^, k = 
1,2,3, ...) be the common solution of the two integrable flows and ^6^), then 

m = ^ k = 1,2,3,..., (67) 

y {Aqi{x,tk),P3{x,tk)f 



satisfy the positive order hierarchy 

In particular, we have the following theorem. 

Theorem 8 Let p{x,t),q{x,t) (p{x,t) = ipi,P2,P3V ,q{x,t) = iqi,q2,q3V) be the 
common solution of the two integrable flows ^6l\) and ^6^), then 

m = —^^^=^^^= (68) 
Y {Aqi{x,t),p3{x,t)f 

is a parametric solution of the 5th-order equation: 

2 2 2 

mt = 4(m"3)^ - 5(m~3)^^^ + {m'^)xxxxx- (69) 
Proof: A direct check is done through substitution of Eqs. ( pTj ) and (pB|). 



6 Peaked stationary solutions 

We know that the equation mt + umx + 3mux = 0, m = u — nixx has peakon solution 
u = e~'^^*L Furthermore, a generalized 6-balanced equation mt + unix + hmux = 
, m = u — nircx , b = constant is found to have this kind of solution ||T^ . 

2 

Now, we study the traveling wave solution of the equation rrit = 4:{m~3)^ — 
5{m~^)xxx + {.'ni~i)xxxxx- Set m~i = v, then this equation becomes 

3 2 

- -V~^Vt=AVx- 5Vxxx + Vxxxxx- (70) 
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Assume this equation has the solution v = f{i), ^ = x — ct, c = constant, then we 
have 

2cr^ = 2f -^r + /'r-^r. (71) 

The right hand side of this equation is quadratically homogeneous and the left hand 
side not. Therefore, we set / = = constant and substitute it into Eq. ([7TD, 

and obtain 

c = 0, (72) 
a^-5a2 + 4 = 0, (73) 

which implies 

a = ±1, a = ±2. (74) 
So, the 5th-order equation ([70| ) has the following stationary solutions: 

1, e-^ e^ e-2^ (75) 
which exactly composes the basis of the solution space of the stationary equation 

2 2 

4:Vx — 5vxxx + Vxxxxx = 0. Therefore, the 5th order PDE mt = 4(m~3 — 5(m~3)^^^ + 
{m~^)xxxxx possesses the stationary solutions 

m{x) = {co + cic-'' + C2e'' + cse-^"" + Cie^"")-^ , (76) 
V Cfc G M, k = 0,1,2,3,4. 

Apparently, e~i^, ei^, e~^^, e^^ satisfy the 5th order PDE nit = 4(m~i)^.— 5(m~i)^.2.^.+ 
{m~^xxxxx- Therefore, /j^(x) = e"tl''l, //(x) = eil'^l, /2^(x) = e"^!"'! , and /2^(a;) = 
g3|a;| g^j^g peaked stationary solutions (see Figure 1) of the 5th-order PDE 

mt = 4(m-i)^ - 5{m-f)xxx + {m~i)xxxxx- 



Comparison with the parametric solution 

All the stationary solutions ([76| ) may be included in the parametric solution ([681). 
For example, the function m{x) = e~2^ is cast in Eq. (^) when we choose dynamical 
variables qi, ps such that the following constraint 

(Agi(x),p3(x)) = e^ (77) 

holds, where gi, are the solutions of the integrable Hamiltonian system (0). But 
the four peaked stationary solutions can not be included in the parametric solution 
( |68D because (Agi(x),p3(x)) is smooth every where. 
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—2/3 —2/3 

Peaked Stationary Solutions for equation mj=4(m )^-5(m )x>,x+(m 




Figure 1: Four peaked stationary solutions 



Comparison with the cusp solution [[22| ] 



Wadati, Ichikawa and Shimizu p2[ proposed the cusp soliton (cuspon) solution for 
the deformed Harry-Dym equation + (1 — rYr^xx = 0. This equation is actually 
equivalent to the Harry-Dym equation qt — 2(q~^^'^)xxx = by the transformation 
1 — r. By the inverse scattering method p|, they obtained the following 



-1/2 



traveling wave solution 



cosh ^, ^ = ax — 4a^t + tanh ^ + A, 



(75 



where a ^ 0,a, A = constants. They called this solution the cusp soliton. It seems 

.-v.^-v.-v. looks like 



5(m 



+ (m~ 



that the present 5th-order PDE mt = A{m~^] 
a higher-order Harry-Dym equation and should have the cusp-kind solution (78). 
But unfortunately, this is not the case for our 5th-order PDE mt 



Aim 



5(m ^i)xxx + ^^)xxxxx- However, our recent study [T^ reveals that the following 
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5th-order PDE 

(15^2 ^ ^^g^ ^ 2)rt + (1 - r)V„ - (79) 
possesses the cuspon r = cosh~^ ^ with ^ — ax — 8a^t + tanh ^ + A, Va 7^ 0, A e R. 
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